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^ Abstract. Hadrontherapy is an emerging technique in cancer therapy that uses 

jy^ beams of charged particles. To meet the improved capability of hadrontherapy 

^~. in matching the dose release with the cancer position, new dose monitoring 

(*^ techniques need to be developed and introduced into clinical use. The 

f^ , measurement of the fluxes of the secondary particles produced by the hadron beam 

I I is of fundamental importance in the design of any dose monitoring device and is 

eagerly needed to tune Monte Carlo simulations. We report the measurements 

^— H done with charged secondary particles produced from the interaction of a 80 

^ MeV/u fully stripped carbon ion beam at the INFN Laboratori Nazionali del 

fvj Sud, Catania, with a Poly-methyl methacrylate target. Charged secondary 

{"^ particles, produced at 90° with respect to the beam axis, have been tracked 

1^^ with a drift chamber, while their energy and time of flight has been measured by 

_^ means of a LYSO scintillator. Secondary protons have been identified exploiting 

• the energy and time of flight information, and their emission region has been 

CO reconstructed backtracking from the drift chamber to the target. Moreover a 

C _ ^ position scan of the target indicates that the reconstructed emission region follows 

^SJ the movement of the expected Bragg peak position. Exploting the reconstruction 

^— H of the emission region, an accuracy on the Bragg peak determination in the 

• • submillimeter range has been obtained. The measured differential production 

1^ rate for protons produced with E^^°'^ > 83 MeV and emitted at 90° with 

kp respect to the beam line is: dNp/{dNcdn){EP[°'^ > 83 MeV, 9 = 90°) = 

rS (2.69 ± 0.08stat ± 0.12sys) X IQ-^sr-i. 
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1. Introduction 

Protons and carbon ion beams are presently used to treat many different solid cancers 
(Jakel et al. 2008, Durante & LoefHer 2010) and several new centers based on hadron 
accelerators are operational or under construction (Amaldi & Kraft 2005, Schardt & 
Elsasser 2010). The main advantage of this technique, in comparison to the standard 
radiotherapy with X-ray beams, is the better localization of the irradiation dose in 
the tumor affected region sparing healthy tissues and possible surrounding organs at 
risk. This feature can be achieved because the heavy charged particles loose most of 
the energy at end of their range, the Bragg peak (hereafter BP), in comparison to the 
exponentially decreasing energy release of the X-ray beam. Up to now most of the 
patients have been treated at centers with proton beams, but routinary use of carbon 
beams has now started. There are also proposals for future use of '^He, ''Li or ^^O 
beams (Brahme 1986). 

New dose monitoring techniques need to be developed and introduced into clinical 
use, to meet the improved capability of hadrontherapy to match the dose release with 
the cancer position. The R&D effort should be then focused to develop novel imaging 
methods to monitor, preferably in real time, the S-dimensional distribution of the 
radiation dose effectively delivered during hadrontherapy. 

This holds true especially for treatments using carbon ion beams since the dose profile 
is very sensitive to anatomical changes and minor patients' positioning uncertainties. 
Conventional methods for the assessment of patients' positioning used in all X-ray 
based radiation therapy, where a non-negligible fraction of the treatment beam is 
transmitted through the patient, cannot be used to pursue this task due to the different 
physics underlying. All the proposed methods exploit the information provided by the 
secondary particles produced by the hadron beam along its path to the tumour, inside 
the patient's body. In particular it has been already shown that the peak of the dose 
released by the hadron beam can be correlated with the emission pattern of the fiux of 
secondary particles created by the beam interaction, namely: i) prompt photons within 
the 1-10 MeV energy range (Min et al. 2006, Testa et al. 2008, Testa et al. 2009) and ii) 
pairs of back-to-back photons produced by the annihilation of positrons coming from 
P+ emitters, mainly "C and ^^O (Pawelke et al. 1997, Parodi et al. 2002, Enghardt 
et al. 2004, Fiedler et al. 2008, Vecchio et al. 2009, Attanasi et al. 2009). 
In this paper we suggest the possibility to correlate the position of the BP in the 
patient with the emission region of charged secondary particles, mainly protons with 
kinetic energy E^inilSO MeV. We report the study of the charged secondary particles 
produced from the irradiation of a Poly-methyl methacrylate (PMMA) with the 80 
MeV/u fully stripped carbon ion beam of the INFN Laboratori Nazionali del Sud 
(LNS). Section [2] is devoted to the description of the experimental setup; the event 
selection and the spectra of the charged secondary particles are presented in Section [3] 
The analysis of the production region of charged secondary particles is described in 
Section HI and the measurement of their differential production rate is reported in 
Section [5] 

2. Experimental setup 

Figure [1] shows the experimental setup. A 4x4x4 cm'^ PMMA target is placed on 
a 80 McV/u, fully stripped ^^C ion beam. The beam rate, ranging from hundreds 
of kHz to ^ 2 MHz, is monitored with a 1.1 mm thick scintillator placed at 17 cm 
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Figure 1. Experimental setup: picture and schematic view. The 
acquisition is triggered by the coincidence of the Start Counter and the 
LYSO detector. 



from the PMMA on the beam line and read-out by two photomultiplier tubes (PMTs) 
Hamamatsu 10583 put in coincidence (Start Counter). 

An array of 4 LYSO crystals, each measuring 1.5x1.5x12 cm^^, is placed at 90° with 
respect to the beam line, at 74 cm from the PMMA center. The scintillation light 
of the crystals is detected with a PMT EMI 9814B triggered in coincidence with the 
Start Counter. 

A 21 cm long drift chamber (Abou-Haidar et al. 2012) is placed at 51 cm from the 
PMMA center, along the flight hne connecting the PMMA to the LYSO crystals. We 
have chosen the configuration at 90° with respect to the beam line to maximize the 
sensitivity to the Bragg peak position along the beam. In the following the coordinate 
system is defined (Figure fTl) with the x-axis along the beam line toward the Start 
Counter, the z-axis along the line connecting the centers of PMMA, drift chamber and 
LYSO crystals detector and the y-axis oriented according to the right-hand rule, and 
origin at the drift chamber center. The PMMA center in the vertical axis was at y = -4 
mm. The drift chamber provides a 2-dimensional reconstruction of the space point by 
alternated horizontal (x-z plane V-view) and vertical (y-z plane U-view) layers of wires. 
Each layer is composed of three 16x10 mm^ rectangular cells for a total of 36 sense 
wires (Figure^. The twelve layers, six on each view, provide tracking redundancy and 
ensure high tracking efficiency and excellent spatial resolution. In order to minimize 
tracking ambiguities, the consecutive layers of each view are staggered by half a cell. 
Custom front-end electronics boards, designed and realized at the INFN Laboratori 
Nazionali di Frascati (LNF) electronics workshop, are embedded in the detector and 
provide single wire signal amplification by a factor of 10. The drift chamber has 
been operated with 1.8 kV sense wire voltage, Ar/C02 (80/20) gas mixture and 30 
mV discriminating threshold for the signals, achieving < 200 /im single cell spatial 
resolution and ~ 96% single cell efficiency (Abou-Haidar et al. 2012). 
The signals from the Start Counter and the LYSO crystals are split and fed into a 
12-bit QDC (Caen V792N) and a 19-bit TDC (Caen V1190B), after discrimination, to 
provide the measurements of both the particles' energy and arrival time. The signals 




Figure 2. Mechanical drawing of the drift chamber. 



from the 36 cells of the drift chamber are fed, after discrimination, into the same TDC 
providing the drift time measurements. The front-end electronics has been read-out 
by a VME system using a MOTOROLA 5100 CPU board. 

The energy and time calibration of the LYSO crystals and the determination of the 
drift chamber space-time relations have been described in (Agodi et al. 2012a, Agodi 
et al. 2012 &) and (Abou-Haidar et al. 2012) respectively. A custom tracking 
algorithm, based on a least squares iterative fitting method, has been also developed 
to reconstruct the direction of the charged secondary particles. A first track 
reconstruction is performed using very clean topologies asking for at least three layers 
with a single fired cell (hit), on both the V- and U- views. Then a hit addition algorithm 
improves the tracking performance by using the information from all the other layers. 

3. Data selection 



The trigger signal is provided by the coincidence of the Start Counter and LYSO 
crystals signals, within 80 ns. Figure [3] shows the distribution of the number of hits 
in the drift chamber ( A'nits ) obtained for events with detected energy in the LYSO 
crystals E'lyso > 1 MeV. Events with A^Hits > 9 are selected for the analysis of the 
charged secondary particles. 

In order to evaluate the setup acceptance and efficiency, and to optimize the particle 
identification analysis a detailed simulation has been developed using the FLUKA 
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Figure 3. Distribution of the number of fired cells in the drift chamber 
A'^Hits for events with detected energy in the LYSO crystals -Blyso > 1 
MeV. 



software release 2011.2 (Fasso' et al. 2003, Ferrari et al. 2005). The detailed geometry 
description with the setup materials (air included) together with the trigger logic, the 
time resolution of the scintillator as well as the experimental space resolution of the 
drift chamber have been considered. The quenching effect in the scintillator has also 
been introduced in the Monte Carlo according to (Wakabayashi et al. 2011). The 
interaction of a sample of 10^ carbon ions with 80 MeV/u, equivalent to 10'^ s of data 
taking at the typical 1 MHz rate of beam, has been simulated. To identify the charged 
particles reconstructed in the drift chamber, we exploit the distribution of the detected 
energy in the LYSO detector E'lyso as a function of Time of Flight (ToF), Figure [4] 
In the data sample (left panel) a fast low-energy component due to electrons is clearly 
visible for ToF values around zero, in the area delimited by the first dashed line. These 
electrons are produced by Compton scattering of the de-excitation photon induced by 
beam interactions in the PMMA material. The central most populated band, delimited 
by the two dashed lines, is made by protons with detected energy within a very wide 
range, originating also the clearly visible saturation of the LYSO crystals QDC for 
EhYso > 24 MeV. The FLUKA simulation (right panel) shows similar populations in 
the (ToF , -Elyso) plane with an additional component of deuterons, above the second 
dashed line, which is not present in data. 

We have then identified as proton a charged secondary particle with ToF and -Elyso 
values inside the area delimited by the two dashed lines in Figure |4j The systematic 
uncertainty on the proton/deuteron identification has been estimated using the data 
events in the deuterons area of the (ToF , E'lyso) plane. Figure Is] shows the 
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Figure 4. Distribution of the detected energy in the LYSO crystals as 
a function of the Time of Fhght: Data (Left) and FLUKA Simulation 
(Right). 



distributions of /3 = - and the corresponding detected kinetic energy -Ekin for the 
identified protons, obtained using the ToF measurement together with the distance 
between LYSO crystals and PMMA. This detected kinetic energy can be related to 
the proton kinetic energy at emission time, E^i^'^, considering the energy loss in the 
PMMA and the quenching effect of the scintillating light for low energy protons. The 
minimum required energy to detect a proton in the LYSO crystals is E^^'^°'^ = 7.0 ±0.5 
MeV, evaluated using the FLUKA simulation, and a proton with an average detected 
kinetic energy E'kin — 60 MeV has been emitted with _E^;™'^83±5 MeV. The uncertainty 
is mainly due to the finite size of both the beam spot 0(1 cm) and profile. 
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Figure 5. Distribution of /3 = ^ (left) and kinetic energy (right) of charged 
secondary particles identified as protons. 



In order to use the secondary protons for monitoring purposes, the crossing of some 
centimeters of patient's tissue has to be considered and therefore the range ii^kin > 
60 MeV of the detected kinetic energy distribution is the most interesting for the 
above-mentioned apphcation. In the following the proton kinetic energy detected in 
the LYSO crystals will be referred to as the kinetic energy. 

4. Production region of charged secondary particles 

Tracks reconstructed in the drift chamber are backward extrapolated to the PMMA 
position, to find the production region of charged secondary particles along the path 
of the carbon ion beam. The PMMA is mounted on a single axis movement stage 
allowing position scans along the x-axis to be performed with a 0.2 mm accuracy 
(Figure fl]). In the configuration with the centers of PMMA, drift chamber and LYSO 
crystals aligned along the z-axis, the PMMA position in the stage reference frame is 
taken as and will be referred to as the reference configuration. 

From each track reconstructed in the drift chamber and backward extrapolated to 
the beam axis we can measure the x and y coordinates of the estimated emission 
point of the charged secondary particle, named xpmma and 2/pmma- The expected 
position of the Bragg peak obtained with the FLUKA simulation (Fasso' et al. 2003) 
is located at (11.0 ± 0.5) mm from the beam entrance face of the PMMA. With the 
setup in the reference configuration, the expected position of the Bragg peak in our 
coordinate system is XBraggl^"^ — (9.0±0.5) mm. Figure^shows the distribution of the 
reconstructed a;pMMA, compared to the expected distribution of the dose deposition in 
the PMMA, both obtained with the setup in the reference configuration. The mean of 
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Figure 6. Expected dose deposition in the PMMA evaluated with FLUKA 
(hatched) compared to the distribution of a;pMMA (solid), the emission 
point of charged secondary particles along the x-axis. The beam entrance 
and exit faces of the PMMA are at a;pMMA = 2 cm and xpmma = -2 cm, 
respectively. 



the gaussian fit to the distribution is xpmma = 17.1 ± 0.2 mm, and consequently 
the separation between the BP and the peak from secondary proton emission is 
AprotonBiagg = 8.1 ± 0.5 mm. Figure [7] shows the distribution of the reconstructed 
a;pMMA and j/pmma for all identified protons (solid line), for protons with i^km > 60 
MeV (hatched) and for protons with i?kin > 100 MeV (grey) . The beam entrance and 
exit faces of the PMMA are at xpmma = 2 cm and xpmma = -2 cm, and j/pmma = 
1.6 cm and ypMMA = -2.4 cm. The xpmma distribution is related to the range of the 
beam while the j/pmma to its transversal profile. Quite remarkably the shape of the 
distribution of the emission point is approximately the same for protons emitted with 
different kinetic energies, e.g. the resolution on xpmma does not depend critically on 
the Skin variable. 
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Figure 7. Distribution of a;pMMA (Left) and 2/pmma (Right) obtained for 
all charged particles identified as protons (black solid line), for protons with 
Skin > 60 MeV (dashed line) and with iJkin > 100 MeV (grey). The beam 
entrance and exit faces of the PMMA are at a;pMMA = 2 cm and a;pMMA 
— -2 cm, and j/pmma ~ 1-6 cm and j/pmma ~ -2.4 cm. 



The existence of a relationship between the expected BP position and the peak of the 
a;pMMA distribution, as a function of the PMMA position, in principle could allow us 
to follow the BP position using the xpmma measurements. To estimate the accuracy of 
this method , a position scan has been performed acquiring several data runs moving 
the PMMA by means of the translation stage. 

For each run with different PMMA position, the production region of the protons 
have been monitored using the mean values of the gaussian fit to xpmma and ypMUA 
distributions, xpmma and j/pmma- Since ypuMA is the coordinate of the proton 
emission point along the vertical axis, and is related to the fixed beam profile in 
the transverse plane, its behaviour as a function of the PMMA position provides an 
estimate of the method's systematic uncertainty. 

Each PMMA position in the stage reference frame can be translated into the expected 
Bragg peak position XBragg for that given PMMA position. Figure [8] shows the results 
obtained for Spmma and ypMMA as a function of a;Bragg, with E'kin > 60 MeV protons. 
A clear linear relationship is observed between Spmma and a^Bragg, indicating that 
the charged secondary particles emission reconstructed with the drift chamber follows 
accurately the BP movement. No dependence of the j/pmma values on the Bragg peak 
position is observed, as expected from a translation of the PMMA along the x-axis 
only. Similar results can be obtained using protons with different i^kin selection, as it 



can be inferred from Figure [7J 
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Figure 8. Reconstructed peak position of the secondary proton emission 
distribution SpmmAjI/pmma as a function of tiie expected Bragg Peak 
position XBragg, with iSkin > 60 MeV. 



To estimate the achievable accuracy on the BP determination several contributions 
need to be considered. We evaluated the difference Api-otonB: 



ragg 



a^PMMA — a^B 



ragg 



for 

all identified protons and for the proton sample with -Ekin > 60 MeV. The AprotonBragg 
root mean square is crAp„tonBra — 0-9 nrm for both samples. This can be explained 
as follows: in the sample with all identified protons the contribution to the total 
uncertainty due to the scattering is partially compensated by the larger statistics with 
respect to the sample with i^kin > 60 MeV. Table [T] reports the number of identified 
protons with E'kin > 60 MeV obtained with the position scan data. 



Table 1. Statistics of identified protons witfi -Bkin > 60 MeV, obtained witfi the 
position scan data. 

a;Bragg (mm) -19 -15 -11 -9-5-3 1 5 9 11 13 

^p?otony 67 77 88 61 92 75 113 154 1223 130 83 



The uncertainty CTExtrapoi due to the backward extrapolation of the track from the drift 
chamber to the beanr line can be estimated from the root mean square of the ypMMA 
values, CTypMMA = cTExtrapoi = 0.5 mm. The latter contributes to the AprotonBragg 
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distribution, together with astago — 0.2 mm from the uncertainty on the PMMA 
positioning. We can then estimate the contribution to the total uncertainty coming 
from the shape of the distribution of the emission point of charged secondary particles 
as: 

^Emission = y^^Ap,ot„„B,a,, " 4xtrapol " ^Stagc ^ O'^mm (1) 

It must be stressed that this value represents only an indication of the precision 
achievable in the BP determination using secondary protons, due to the target 
thickness and homogeneity in the present setup, with respect to a possible clinical 
application. 

5. Flux of charged secondary particles 

The flux of the secondary protons emitted from the beam interaction with the PMMA 
has been measured at 90° with respect to the beam direction and in the geometrical 
acceptance of the triggering LYSO crystals, configuration maximizing the sensitivity 
to the Bragg peak position. The surface of the LYSO is 3x3 cm^, corresponding to 
a solid angle ilLYSO = 1-3 x 10~^ sr at a distance of 74 cm. The proton's kinetic 
energy spectrum measured with data has been inserted in the FLUKA simulation to 
evaluate the detection efficiency in the LYSO crystals for protons with E'lyso > 1 
MeV: elyso = (98.5 ± 1.5)%, with the uncertainty mainly due to the Monte Carlo 
statistics. To properly evaluate the rate of charged secondary particles reaching the 
LYSO crystals, the number of carbon ions reaching the PMMA target (A^c) has 
been computed according to (Agodi et al. 2012a): counting the number of signals 
in the Start Counter (A^sc) within randomly-triggered time-windows of T^, = 2 |j,s, 
corrected for the Start Counter efficiency esc = (96 ±1)%, and the acquisition 
dead time. The number of emitted secondary protons A'p has been measured with 
the a;pMMA distribution counts, corrected for esci clysOi the tracking efficiency 
CTrack = (98 ± 1)% (Abou-Haidar et al. 2012) and the acquisition dead time. 
The double differential production rate of secondary protons emitted at 90° with 
respect to the beam line is estimated as: 

dNcdn^'-^'' >-Nc f^LYSo' ^ ^ ' 

Figure [9] shows the double differential production rate of secondary protons, emitted 
at 90° with respect to the beam line, as a function of the rate of the carbon ions Rq 
reaching the PMMA: all identified protons and protons with E^in > 60 MeV. 
Expressing these results in terms of the secondary proton's kinetic energy at emission 
^kin°'^' we obtain: 

"^^P iE^[°'' > 7 MeV, 9 = 90°) = (9.56 ± 0.18stat ± 0.40sy.) x lO-^sr-^ (3) 



dNcdn 



-^^«i;°'^ > 83 MeV, - 90°) - (2.69 ± 0.08,tat ± 0.12^^^) x lO-^sr-^ (4) 

aiVcui' 

with the systematic contribution mainly due to proton identification and to the 

uncertainty on the production kinetic energy related to the beam's transversal profile 

uncertainty. 

The same experimental setup described in Section |2] has been used to measure the 

differential production rate for prompt photons, with energy i?LYSO > 2 MeV and 

emitted at 90° with respect to the beam line: dN^/{dNcd^){Ei^Yso > 2 MeV, 6* = 

90°) = (2.92 ± 0.19) X lO-^sr-i (Agodi et al. 2012a). 



11 



v_ 0.0012 


— 








-!2- 


- - 










^o 0.001 




A 




.. , , 1 


i 


k 


i 








3 




J fc 


T-^ 


' - 11 


1 0.0008 


— 




t 




1 


0.0006 


— 








0.0004 


- 




J 


L 




= *A * 


Tj 


^ ^ 


L ^ 




^ • 








' T 




_ T 




9 


i I -^ 


0.0002 


-,,,!,,,! 


#* 


* 

,!,,,! 


1 1 1 1 1 1 1 1 1 1 



0.8 



1.2 



1.4 



1.6 



1.8 



Re (MHz) 



Figure 9. Double differential production rate for secondary particles 
emitted at 90° with respect to the beam line, as a function of the rate 
of the carbon ions Re reaching the PMMA target: all identified protons 
(triangles) and protons with Skin > 60 MeV (circles). 



6. Discussion and conclusions 



We reported the study of secondary charged particles produced by the interaction of 
80 MeV/u fully stripped carbon ion beam of INFN-LNS laboratory in Catania with a 
PMMA target. Protons have been identified exploiting the energy and time of flight 
measured with a plastic scintillator together with LYSO crystals, and their direction 
has been reconstructed with a drift chamber. A detailed simulation of the setup based 
on the FLUKA package has been done to evaluate its acceptance and efficiency, and 
to optimize secondary particle's identification. 

It has been shown that the backtracking of secondary protons allows their emission 
region in the target to be reconstructed. Moreover the existence of a correlation 
between the reconstructed production region of secondary protons and the Bragg peak 
position has been observed, performing a position scan of the target. The achievable 
accuracy on the Bragg peak determination exploting this procedure has been estimated 
to be in the submillimeter range, using the described setup and selecting secondary 
protons with kinetic energy at emission E^^°'^ > 83 MeV. 

The obtained accuracy on the position of the released dose should be regarded as 
an indication of the achievable accuracy for possible applications of this technique to 
monitor the BP position in hadrontherapy treatment. In fact in clinical application the 
secondary particles should cross a larger amount of material (patient tissue) resulting 
in an increased multiple scattering contribution worsening the BP resolution by, at 
most, a factor 2-3. On the other hand an optimized device allowing a closer positioning 
to the patient could greatly improve the collected statistics of protons produced with 
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E^[°'^ > 80 McV, reducing multiple scattering effects. Furthermore the intrinsic good 
tracking resolution and high detection efficiency easily achievable in charged particles 
detectors, make this monitoring option worthwhile of further investigations. 
The measured differential production rate for protons with -E'^;™'^ > 83 MeV and 
emitted at 90° with respect to the beam line is: dNp/{dNcdn){E^{°'^ > 83MeV,e ^ 
90°) == (2.69 ± 0.08stat ± 0.12sys) x lO-*sr-\ 
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